Topological Properties of Phase Singularities in Wave Fields 
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Phase singularities as topological objects of wave fields appear in a variety of physical, chemical, 
and biological scenarios. In this paper, by making use of the fli-mapping topological current theory, 
we study the topological properties of the phase singularities in two and three dimensional space in 
details. The topological inner structure of the phase singularities are obtained, and the topological 
charge of the phase singularities are expressed by the topological numbers: Hopf indices and Brouwer 
degrees. Furthermore, the topological invariant of the closed and knotted phase singularities in three 
dimensional space are also discussed in details. 
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I. INTRODUCTION 

In physics, the solutions of wave equations in 
two and three dimensional space often possess phase 
singularities [J 0]. At the point of phase singularity, the 
phase of the wave is undefined and wave intensity van- 
ishes. This phase singularity is a common property to 
all waves, it exists in many area of physical, chemical, 
and biological scenarios, such as the quantized vortices 
in superfluid or superconductor systems 0, Q, the vor- 
tices phase singularities in Bose-Einstein condensates @, 
the streamlines or singularities in quantum mechanical 
wavefunctionsjaja, th e op tical vortices in optical wave 
systems l9l. Hoi [ill [l2l fl3| . the vortex filaments in chem- 
ical reaction and molecular diffusion [IH, [l6l 171 . and 
the phase singularities in biological systems 17 , lg] . In 
recent years, the phase singularities have drawn great 
interest because it is of importance for understanding 
fundamental wave physics and have many important 
applications, and a great deal of works on the phase 
singularities in wave fields have been done by many 

physicistsfE H, HE HI - 

Phase singularities are generically points in two dimen- 
sional space and form a network of lines in three di- 
mensional space. Around the phase singularity points 
in two dimensions or phase singularity lines in three 
dimensions, the phase changes by 2tt times an inte- 
ger. This integer is the so-called topological charge (the 
strength of the singularity) of the phase singularity ob- 
jects (points or lines). In three dimensional space, a im- 
portant case is that the phase singularity lines are closed 
and knotted curves, these knot-like configurations exist 
in a variety of physical scenarios, including Bose-Einstein 
condensatio ns 1141 . c hemical reaction and molecular diffu- 
sion system s Il5l. 1 16l [l~7| . optical wave systems [l9l,l20j and 
field theory [2ll |22| | . The phase singularities (points, lines, 
or knotted) are topological objects, and have a very rich 
topological properties, so it is necessary to use the topo- 
logical method to study the topology of these objects. 



In this paper, we have given a generic feature of the 
topological properties of the phase singularities in wave 
fields. By making use of the 0-mapping topological 
current theory [23l|, we study the topological current of 
the phase singularity objects (points in two dimensions 
and lines in three dimensions), and the topological inner 
structure of these topological objects are obtained. For 
the case that the phase singularity lines in three dimen- 
sional space are closed and knotted curves, we discussed 
the topological invariant of these knotted family in de- 
tails. This paper is arranged as follows. In Sec. II, we 
construct a topological current density of the phase sin- 
gularity points in two dimensions, this topological cur- 
rent density don't vanish only when the phase singular- 
ity points exist, the topological charge of phase singu- 
larity points are expressed by the topological quantum 
numbers, the Hopf indices and Brouwer degrees of the <f>- 
mapping. In Sec. Ill is the topology of the phase singu- 
larity lines in three dimensions. In Sec. IV, we introduce 
a important topological invariant to describe the phase 
singularity lines when they are linked and knotted, it is 
just the sum of all the self-linking and all the linking 
numbers of the knot family. In Sec.V is our concluding 
remarks. 



II. TOPOLOGICAL STRUCTURE OF PHASE 
SINGULARITIES IN TWO DIMENSIONS 

Let us study a complex scalar field ip(x,t), the wave 
fields function ip(x,t) is maps from space (of either 
two or three dimensions) to the complex numbers, so 
ijj : R 2 , R 3 — ► C. The complex scalar field ip(x, t) is time 
dependent as well as space dependent, and it can be writ- 
ten as 



i/>(2,t) = HVII e l9 =4> 1 {x,t)+i<t> 2 {x,t), 



(1) 
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where (f) 1 (x, t) and 4> 2 (x, t) can be regarded as complex 
representation of a two-dimensional vector field (f> — 
{(f) 1 , (f) 2 ) over two or three dimensional space, ||^|| = 
^ip*ip is the modulus of yj, and 9 is the phase factor. 
In this section, we will restrict our attention to the phase 
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singularities in two dimensional space, and labeled points 
x = (a; 1 , a; 2 ) in cartesian coordinates. The phase singu- 
larities in three dimensional space will be discussed in 
Sec.III and Sec.IV. 

It is known that the current density J associated with 
4>{x,t) is defined as 



J = 7m(V>*VV0 = U\\ 2 v, 



(2) 



the v is the velocity field. From the expressions in Eq. |T|) 
and Eq.©, the velocity field v can be rewritten as 



1 1 

2* PI 



r O*vv>- W>» =V0, 



(3) 



it becomes the gradient of the phase factor 8. The vor- 
ticity field uj of the velocity field v is defined as 



LO = — V X V. 

2vr 



(4) 



Obviously, in two dimensional space the uj only has z- 
component, i.e., uj = uj ■ e z . Form Eq.©,we directly 
obtain a trivial curl-free result: uj — ^-V x V9 = 0. But 
in topology because of the existence of phase singularities 
in the wave fields ip, the vorticity uj does not vanish]^]. 
So in the following discussions, we will study that what 
the exact expression for uj is in topology. 

Introducing the unit vector n a = a /||</>||(a = 
l,2;n a n a — 1), and defining a topological current in 
(2+1) dimensions space-time, 



Ztt 



1 = 0,1,2. (5) 



Obviously, the topological current is identically con- 
served, 



dd l = 0. 



(6) 



By making use of the definition of the unit vector n a , 
one can reexpressed the velocity field v as 



r* an b 

v = e ab n Vn , 



and the vorticity uj is 

.n 1 



2tt 



e^ k e ab d ? n a d k 



(J) 



J,* = 1,2, (8) 



it is just the time component of the topological cur- 
rent f . According to the </>-mapping topological current 
theory [23j], one can prove that 



where 



X 



D*{£) - le ijk e ab d j( f> a d k <f> b 



(9) 



(10) 



is the Jacobian vector, in which the time component 
D°(4>/x) is the usual two dimensional Jacobian D((f>/x). 



The expression of the topological current in Eq. ([9]) shows 
that the topological current j l does not vanish only at 
the zero points of <fi, i.e., the phase singularities of the 
wave fields ip. According to the Eq.© and Eq.©, the 
vorticity can be expressed 



u J =f=5\4>)D^), 
x 



(11) 



obviously, uj also does not vanish only when the phase 
singularities of the wave fields ip exist. The topological 
current j 1 and vorticity uj paly an important role in the 
topological property of the phase singularities, so it is 
necessary to study the zero points of <p to determine the 
non-zero solutions of j % . 

According to the implicit function theory|24| , while the 
regular condition 

is satisfied, the general solutions of 

^(x^x 2 ,^ = 0, <f> 2 (x 1 ,a?,t)=0, (12) 

can be expressed as 

x l =xl(s,t), x 2 = x 2 k {s,t), k = l,2,...,N, (13) 

which represent N zero points zu (k=l, 2, . . . , N) where 
ip(f,t) = in two-dimensional space. These zero points 
solutions are just the so-called phase singularity points 
of the wave fields ip, they represent the zeroes of the real 
part and imaginary part of the complex scalar wave fields: 
Re(tp) = (j) 1 = and Im(ip) = 2 = 0. 

In the theory of i5 function of 4>(r), one can prove that 
in 2-dimensional space [25| 



5 2 {dj) 



N 

E 



i \DWx)g h 



■5 2 (x- z k ), 



(14) 



where the positive integer (3k is called the Hopf index of 
map x — > </i[l^. [l3|. The meaning of (3k is that when the 
point z*k covers the neighborhood of zero Zk once, the 
vector field (j) covers the corresponding region (3k times. 
By substituting Eq. (fl"4|) into Eq.®, one can obtain that 



dx i 



N 

f = Y / WkS 2 (x~Zk)^ F , 
fe=i 



(15) 



where Wk = (3kT\k is the winding number of <j) around zero 
point Zk, and r\k = sgn(D((f>/x)) = ±1 is the Brouwer 
degrees of </>-mapping. The sign of Brouwer degrees are 
very important, for the case of vortices phase singulari- 
ties, the rjk = +1 corresponds to the vortex, and rjk = — 1 
corresponds to the antivortex. 

The topological charge (the strength of phase singular- 
ities) of the phase singularities point Zk is defined by the 
Gauss map n : <9£fc — > 5 ,1 p6l|: 



Wk 



1 

27 



n*(e ab n a dn b 



(16) 



<9£ fc 
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where n* is the pullback of Gauss map n, dT,k is the 
boundary of a neighborhood of point z). and Sfc n 
S m = for E m is the neighborhood of another arbitrary 
zero point z m . In topology it means that, when the point 
z\ covers d£/c once, the unit vector n will cover S , or 
(/> covers the corresponding region Wk times, which is a 
topological invariant. By using the Stokes' theorem and 
in term of Eq.©, one can obtain that 

W k = ±- / e ab e^d 3 n a d k n b d 2 x 

tod 2 x = [3kr]k- (17) 

It is clear that the topological charges densities of phase 
singularity points is just the non- vanishing vorticity to of 
the velocity field v. Using the notion j = (io 7 j), Eq.© 
can be rewrite as 



~dt 



V • j = 0, 



(18) 



this expression tells us that the topological charges of the 
phase singularity points are conserved. This is the only 
topological property of the complex wave fields. 

In this section, we obtain the topological inner struc- 
ture of the phase singularity points in two dimensional 
space, the solutions of the Eq. (fT3"l) possess N isolated 
phase singularity points of which the fcth point possesses 
topological charge W k = PkVk- 



III. TOPOLOGICAL CURRENT OF PHASE 
SINGULARITIES IN THREE DIMENSIONS 

In this section, we will study the phase singularities in 
three dimensions space, some results which discussed in 
the above section in two dimensions space are also useful 
here. Differently with in two dimensions, the point is 
labeled x — (x 1 , x 2 , x 3 ) and the vorticity Co have its all 
three components. In three dimensional space, the Co can 
be reexpressed as 



1 



w 1 = —e t3 *e ab d j n a d k n 1 ', 



(19) 



it is clearly that the vorticity Co is just the topological 
current of phase singularities in three dimensions space. 

Using the </>-mapping theorv(23j|. the topological cur- 
rent Co is rewritten as 



CO 



5 2 {<i>)D\-\ 
x 



(20) 



we can see from this expression that the vorticity Co is 
non-vanishing only if <f> = 0, i.e., the existence of the 
phase singularities, so it is necessary to study these zero 
solutions of 4>. In three dimensions space, these solutions 
are some isolate zero lines, which are the so-called phase 
singularity lines in three dimensions space. 



Under the regular condition 

D l (<P/x) 0, 
the general solutions of 

4> 1 (x l ,x 2 ,x 3 ,t) = 0, <j?(x 1 ,x 2 ,x 3 ,t)=Q (21) 

can be expressed as 

x 1 = xl(a,t), x 2 = x 2 (s,t), x 3 = x 3 k ( S ,t), (22) 

which represent the world surfaces of N moving iso- 
lated singular strings Lk with string parameter s (k — 
1, 2, • • • ,N). These singular strings solutions are just the 
phase singularities solutions in three dimensions space. 

In S- function theory (25j. one can obtain in three di- 
mensions space 



k=l 



(23) 



where 



~ 2 £ €mn dm du k ' 

and Eft is the fcth planar element transverse to Lk with 
local coordinates (u^u 2 ). The (3k is the Hopf index of 
<f> mapping, which means that when x covers the neigh- 
borhood of the zero point c?fc(s) once, the vector field 
4> covers the corresponding region in 4> space [3k times. 
Meanwhile the direction vector of Lk is given by 



dx l 



D^/u) 1 



(24) 



Then from Eq. (|23|) and Eq. (|24|) one can obtain the inner 
structure of to" 1 : 



N 

E 



dr l 

W k I ^S 3 (x-x k (s))ds, (25) 



where Wk = (3kVk is the winding number of <f> around Lk, 
with rjk = sgnD((f)/u)\s k — ±1 being the Brouwer degree 
of <fi mapping. The sign of Brouwer degrees are very im- 
portant, for the case of vortices phase singularities, the 
r]k = +1 corresponds to the vortex, and rjk = —1 corre- 
sponds to the antivortex. Hence the topological charge 
of the phase singularity line Lk is [21] 



io % d<Ji = Wk- 



(26) 



The results in this section show us the topological inner 
structure of the topological current density lo 1 . The topo- 
logical charge of the fcth phase singularity line in three 
dimensions can be expressed by the topological numbers: 
Qk = Wk = /3 k r]k- 
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IV. TOPOLOGICAL ASPECT ON KNOTTED 
PHASE SINGULARITY LINES 

An important case is that the phase singularity lines in 
three dimensions space from closed and knotted curves. 
Topology has play a very important role in understanding 
these knot configurations, so it is necessary to study the 
topology in the knotted phase singularity lines. In order 
to do that, we define a the helicity integral (27) 



H = / r • V .- nt'.v. 



(27) 



this is an important topological knot invariant and it 
measures the linking of the phase singularity lines. From 
the Eq.Q, the helicity integral can be changed as 



H = — I v ■ ujd 3 x. 



1 

2^ 

Substituting Eq. ([25]) into Eq. ([28]). one can obtain 



(28) 



H 



i N r 



■ dx. 



(29) 



when these phase singularities are closed and knotted 
lines, i.e., a family of knots £k(k = 1,2,..., N), Eq. |29|) 
becomes 



1 N 
2tt f-^ 



fc=i 



v ■ dx. 



(30) 



It is well known that many important topological num- 
bers are related to a knot family such as the self-linking 
number and Gauss linking number. In order to discuss 
these topological numbers of knotted phase singularity 
lines, we define Gauss mapping: 



m : S 1 x S 1 S 2 , 
where m is a unit vector 



m(x, y) 



\y- x\ 



(31) 



(32) 



where x and y are two points, respectively, on the knots 
and (in particular, when x and y are the same point 
on the same knot £ , n is just the unit tangent vector T 
of £ at x ). Therefore, when a; and y , respectively, cover 
the closed curves and £j once, n becomes the section 
of sphere bundle 5 2 . So, on this S 2 we can define the 
two-dimensional unit vector e = e(x, y). e, fa are normal 
to each other, i.e. , 

e l • G2 = e i • fa = e*2 • fa = 0, 

e% ■ e% = €2 ■ &2 = m • m = 1. (33) 

In fact, the velocity field u can be decomposed in terms of 
this two-dimensional unit vector e: = e a be a die b — (9;$, 



where 9 is a phase factor Q. Since one can see from the 
expression ui — -^V x v that the (8$) term does not 
contribute to the integral H, Vi can in fact be expressed 
as 

Vi = e ab e a d ie b . (34) 
Substituting it into Eq.(14), one can obtain 

n . 

H = —Y,W k j e ab e a (x , y)d ie b (x, y)dx\ (35) 

Noticing the symmetry between the points x and y in 
Eq. (f32"|) , Eq. ([35|) should be reexpressed as 



H = T- E W kWl £ i e ab d i e a d j e h dx i A dyK (36) 



k,l=l 



In this expression there are three cases: (1) £fc and are 
two different phase singularities (£& ^ £z), and x and y 
are therefore two different points (x ^ y); (2) £/. and 
are the same phase singularities = £/), but x and y 
are two different points (x 7^ y); (3) and are the 
same phase singularities = and x and y are the 
same points (x = y). Thus, Eq. (|36[) can be written as 
three terms: 



H = 



N 1 



■wt 



N 



e a bdie a dje b dx l A dy- 7 



+ 7^ Vwf / e afc e a d 4 e fc dx 4 

iV 

+ E J-^fcWj / j e ab d i e a d j e b dx i Ady J .(37) 

By making use of the relation e a bdie a dje b — ^fa ■ (diffi x 
djfA)^, the Eq. ([37|) is just 



fa*(dS) 



4tt fc 

1 ^ /• 

271 t{ J** 

E ^ WkWl j c j ™*(^)> (38) 



k,l=l(k^l) 



where m* (dS) = m • (dim x djfn)dx l A dy- 7 (x ^ y) denotes 
the pullback of the S 2 surface element. 

In the following we will investigate the three terms in 
the Eq.® in detail. Firstly, the first term of Eq.([3"8) is 
just related to the writhing number [29]| Wr(£k) of 



Wrfa) = ^~<f <f m*(dS) 



(39) 
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For the second term, one can prove that it is related to 
the twisting number Tw(^ k ) of 



I e ab e a d i e b dx i = 



(40) 



where T is the unit tangent vector of knot at x (m = T 
when x = y) and V is defined as e° = e ab V b (V _L f , e = 
fx V). In terms of the White formula^ 



SL(£ k ) = Wr(£ k )+Tw(S k ), 



(41) 



we see that the first and the second terms of Eq. p8|) just 
compose the self-linking numbers of knots. 

Secondly, for the third term, one can prove that 



i- 6 6 m*(dS) 



47T 



1 

47r' 



.i)k 



dx l 



dy 



i (x k - y k ) 
\x-y\\ 3 



(42) 



where Lfc(£fc,£/) is the Gauss linking number between 
and ^(H. Therefore, from Eqs.jH, (|40]), gl|) and (|42|. 
we obtain the important result: 

AT AT 

^ = E W/ ^ L ^ fe ) + E W*W,Lfc(&,£,). (43) 

fc=l k,l=l(k^tl) 

This precise expression just reveals the relationship be- 
tween 7? and the self-linking and the linking numbers of 
the phase singularity knots family 29j. Since the self- 
linking and the linking numbers are both the invari- 
ant characteristic numbers of the phase singularity knots 
family in topology, H is an important topological invari- 
ant required to describe the linked phase singularities in 
wave fields. 



CONCLUSION 



In the present study, by making use of the ^-mapping 
topological current theory, the topology of the phase sin- 
gularity in wave fields is studied. Firstly, we obtain the 
inner structure of the phase singularities in two and three 
dimensional space. The phase singularity objects have 
been found at the every zero point of the wave func- 
tion tp under the condition that the Jacobian determi- 
nate D l {(j)/x) 7^ 0. One shows that the topological charge 
(strength of the phase singularity) of the phase singular- 
ity objects are determined by the topological numbers: 
Hopf indices and Brouwer degrees. Secondly, we have 
studied the topological invariant of the knotted phase sin- 
gularity lines in three dimensional space in details. This 
topological invariant can be expressed as the sum of all 
the self- linking and all the linking numbers of the knotted 
phase singularity lines. 

Finally, it should be pointed out that in the present pa- 
per, when we discussed the topological properties of the 
phase singularities, the condition D l (<fi/x) ^ must be 
satisfied. Now the question is coming, when this condi- 
tion fails, what will happen about the phase singularity 
objects? The answer is related to the evolution of the 
phase singularity objects 31]. The topological object will 
generate, annihilate, split, or merge when the condition 
fails , and these dynamics properties of the phase singu- 
larity objects will be discussed in further work. 



Acknowledgments 



This work was supported by the National Natural Sci- 
ence Foundation of China. 



[1] L. M. Pismen, Vortices in Nolinear fields, (Clarendon 
Press, Oxford, 1999). 

[2] J. Adachi, G. Ishikawa, arxiv: |math/0608721| 

[3] R. J. Donnelly, Quantized Vortices in Helium II (Cam- 
bridge University Press, Cambridge, England, 1991). 

[4] Y. S. Duan, X. Liu and P. M. Zhang, J. Phys: Condens. 
Mather 14, 7941 (2002). 

[5] S. Inouye, S. Gupta, T. Rosenband, A. P. Chikkatur, 
A. Gorlitz, T. L. Gustavson, A. E. Leanhardt, D. E. 
Pritchard, and W. Ketterle, Phys. Rev. Lett. 87, 080402 
(2001). 

[6] J. O. Hirschfelder, A. C. Christoph, and W. E. Palke, J. 

Chem. Phys. 61, 5456 (1974). 
[7] J. Riess, Phys. Rev. D 2, 647 (1970); J. Phys. A: Math. 

Gen. 20, 5179 (1987). 
[8] A. S. Desyatnikov, Y. S. Kivshar, and L. Torner, Optical 



Vortices and Vortex Solitons, Progress in Optics Vol. 47 
(North-Holland, Amsterdam, 2005). 
[9] J. F. Nye, Natural focusing and fine struture of light 
(Institute of Physics Publishing, 1999) . 
[10] M. V. Berry and M. R. Dennis, Proc. R. Soc. A 456, 2059 
(2000). 

[11] M. V. Berry and J. F. Nye, Proc. R. Soc. A 336, 165 
(1974). 

[12] Y. S. Duan, G. Jia and H. Zhang, J. Phys. A: Math. Gen. 

32, 4943 (1999). 
[13] Y. S. Duan, H. Zhang and G. Jia, Phys. Letts. A 253, 57 

(1999). 

[14] E. Babaev, L. D. Faddeev, and A. J. Niemi, Phys. Rev. 
B 65, 100512 (2002); E. Babaev, Phys. Rev. Lett. 88, 
177002 (2002). 

[15] A. T. Winfree and S. H. Strogatz, Physica D 8, 35 (1983); 



6 



Physica D 9, 65 (1983); Physica D 9, 333 (1983); Physica 
D 13, 221 (1984); Physica D 17, 109 (1985). 

[16] A. T. Winfree, Physica D 84, 126 (1995). 

[17] A. T. Winfree, Nature 371, 233 (1994); Science 266, 1003 
(1994). 

[18] A. T. Winfree, The Geometry of Biological Time 

(Springer, 1980). 
[19] M. V. Berry and M. R. Dennis, Proc. R. Soc. A 457, 2251 

(2001). 

[20] J. Leach, M. R. Dennis, J. Courtial and M. J. Padgett, 
Nature 432, 165 (2004), New. J. Phys. 7, 55 (2005). 

[21] Y. S. Duan, X. Liu, and L. B. Fu, Phys. Rev. D 67, 
085022 (2003). 

[22] F. Faddeev and A. J. Niemi, Nature 387, 58 (1997). 

[23] Y. S. Duan, SLAC-PUB-3301, 1984; Y. S. Duan, H. 
Zhang, G. Jia, Phys. Letts. A 253, 57 (1999). 



[24] E. Goursat, A Course in Mathematical Analysis, trans- 
lated by E. R. Hedrick (Dover, New York, 1904), Vol.1. 

[25] J. A. Schouten, Tensor Analysis for Physicists (Claren- 
don Press, Oxford, 1951). 

[26] Y. S. Duan, S. Li and G. H. Yang, Nucl. Phys. B 514, 
705 (1998). 

[27] H. Aref and I. Zawadzki, Nature 354, 50 (1991). 

[28] N. D. Mermin, T. L. Ho, Phys. Rev. Lett. 36,594 (1976). 

[29] E. Witten, Commun. Math. Phys. 121, 351 (1989); A. M. 

Polyakov, Mod. Phys. Lett. A 3, 325 (1988). 
[30] D. Rolfsen, Knots and Links (Publish or Perish, Berkeley, 

CA, 1976). 

[31] L. B. Fu, Y. S. Duan, and H. Zhang, Phys. Rev. D 61, 
045004 (2000). 



